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Abstract 

This report presents the use of a time-marching three- 
dimensional (3-D) compressible Navier-Stokes equation 
numerical solver with a one-equation turbulence model to 
simulate the flow fields developed along concave wall 
surfaces without and with a downstream extension flat wall 
surface. The 3-D Navier-Stokes numerical solver came from 
the NASA Glenn-HT code. The one-equation turbulence 
model was derived from the Spalart and Allmaras model. 
The computational approach was first calibrated with the 
computations of the velocity and Reynolds shear stress 
profiles of a steady flat plate boundary layer flow. The 
computational approach was then used to simulate 
developing boundary layer flows along concave wall 
surfaces without and with a downstream extension wall. The 
author investigated the computational results of surface 
friction factors, near surface velocity components, near wall 
temperatures, and a turbulent shear stress component in 
terms of turbulence modeling, computational mesh 
configurations, inlet turbulence level, and time iteration step. 
The computational results were compared with existing 
measurements of skin friction factors, velocity components, 
and shear stresses of the developing boundary layer flows. 
With a fine computational mesh and a one-equation model, 
the computational approach could predict accurately the skin 
friction factors, near surface velocity and temperature, and 
shear stress within the flows. The computed velocity 
components and shear stresses also showed the vortices 
effect on the velocity variations over a concave wall. The 
computed eddy viscosities at the near wall locations were 
also compared with the results from a two equation 
turbulence modeling technique. The inlet turbulence length 
scale was found to have little effect on the eddy viscosities at 
locations near the concave wall surface. The eddy viscosities, 
from the one-equation and two-equation modeling, were 
comparable at most stream-wise stations. The present one- 
equation turbulence model is an effective approach for 
turbulence modeling in the near solid wall surface region of 
flow over a concave wall. 

Introduction 

The Gas Turbine Branch of NASA Glenn Research Center 
has developed the NASA Glenn-HT CFD Code for the 
prediction of gas turbine blade surface heat transfer. This 


Code used time-stepping Reynolds Averaged Navier-Stokes 
computational method (ref. 1) with a two-equation k - co 
turbulence model (ref. 2) to simulate a turbine blade flow. 
Spalart and Allmaras (ref. 3) assembled an incompressible 
flow one-equation eddy viscosity model based on the 
turbulence modeling of two-dimensional (2-D) free shear 
layer, near wall finite/high Reynolds flows, and transition 
and tripping of laminar flow. With these assembling 
procedures, the present author believes that the one-equation 
eddy viscosity model has the potential for turbulence 
modeling of the Navier-Stokes numerical simulation of near 
surface flow along a turbine blade surface. Before the model 
could be used for turbulence modeling of a complex turbine 
blade flow, the turbulence model should be verified for its 
effectiveness to aid the Navier-Stokes computations of 
simple 3-D boundary layer flows. The turbulence intensity 
and length scale at the inlet of turbine blade and blade 
surface curvature play an important role in the flow 
development along the blade wall surface. The author has 
previously studied (refs. 4 to 6) the applications of the one- 
equation turbulence model and a 3-D Navier-Stokes 
numerical solver (ref. 1) for the simulations of boundary 
layer flow developments related to turbine blade flow field. 

Based on the Spalart and Allmaras eddy viscosity model 
(ref. 3), this author derived a 3-D compressible one-equation 
eddy viscosity equation (ref. 4). The equation was written in 
a form similar to the compressible Navier-Stokes equation. 
The eddy viscosity equation was implemented into the 
computational procedures of a version of the NASA Glenn- 
HT code. Numerical computations were performed (ref. 4) to 
calculate the mean velocity of flat plate boundary layer flows 
with and without flow injection through discrete holes on the 
plate surface. The computational approach predicted the 
occurrence of a fully developed turbulent boundary layer 
velocity profile within the flow over a flat plate surface. The 
approach could also simulate successfully the effect of the 
injection fluid turbulence level on the mean velocity profile 
near an injection hole. Subsequently, this author used the 
computational approach to analyze the development of 
boundary layer flows along concave wall surfaces without 
(ref. 5) and with (ref. 6) a downstream recovery flat wall 
surface. With a global time step (ref. 1) used in these 
previous computational flow simulations, the computed 
friction factors, velocity, and the static temperatures at the 
near wall surface locations oscillated periodically in terms of 
the time iteration step and physical locations along the span- 
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wise direction. These results agreed well with similar results 
from existing experiments (refs. 7 to 9). The numerical 
approach was also found capable of capturing the vortices 
effect on the velocity variation and flow transition along a 
concave wall surface (ref. 7). 

This paper describes the author’s recent work to 
investigate further his previous observations of the results in 
references 4 to 6. The work explored the computational 
predictions, using the one-equation turbulence model (ref. 4), 
of the Reynolds shear stress component, the near wall surface 
temperature, and the occurrence of flow transition as the flow 
developed along a concave wall surface. The near wall eddy 
viscosities from the one-equation turbulence model were also 
compared with the results from a low Reynolds number two- 
equation turbulence modeling technique (ref. 10). First, the 
computed velocity and eddy viscosity profiles (ref. 4) were 
used to verify the ability to calculate the Reynolds shear 
stresses of a fully developed turbulent flat plate boundary 
layer flow. Then, a fine computational mesh was used to 
compute the skin friction factor, the Reynolds shear stress, 
the near wall surface temperature, and the mean flow 
velocity components within the developing flows over 
concave wall surfaces without and with a downstream 
extension wall surface. The computed skin friction factors, 
velocity components, near wall surface temperatures, and a 
component of Reynolds shear stress were compared with 
their existing results in the previous work (refs. 5 and 6) and 
in experiments (refs. 7 to 9). The comparisons showed that 
the accuracy of the numerical computations depended on the 
computational mesh configuration. With a fine mesh, the 
present computational approach could predict accurately the 
variations of velocity, Reynolds shear stress, near wall 
surface temperature, and the occurrence of flow transition. 
The present one-equation turbulence modeling did not 
consider the effect of the inlet turbulence length scale on the 
numerical computations. Alternately, the author used a low 
Reynolds number k-a> turbulence model (ref. 10) to compute 
the eddy viscosities within the flow over the concave wall 
surface. The author compared the eddy viscosities obtained 
from both the one-equation and two-equation turbulence 
modeling. The inlet turbulence length scale was found to 
have little effect on the turbulence modeling properties. The 
near wall surface eddy viscosities from the one-equation 
turbulence model were comparable to the eddy viscosities 
from the two equation turbulence modeling. Thus, the 
present one-equation turbulence model is an effective 
approach to model the turbulence eddy viscosities at near 
solid wall surface locations within a flow developed along a 
concave wall surface. 

Nomenclature 

Cf friction factor 

C p specific heat 

dV s shear layer velocity defect 


E, total energy 

h turning duct height, 0.375 ft 

K heat conductivity 

k turbulence kinetic energy 

L reference length, 0.833 ft or h 

N coordinate normal to wall surface 

m flow rate 

n number of time iteration steps 

p pressure 

Pr Prandtl number 

Q w wall surface heat transfer rate 

R universal gas constant 

Re Reynolds number/unit length, p Ujl p 

Re/. Reynolds number based on L 

Re v Reynolds number based on x 

S Stream-wise coordinate, slh 

T temperature 

Uf free stream velocity 

u, v, w velocity components in x, y, z directions 

V„, V, velocity components in N, S directions 

v working variable for eddy viscosity 

v, eddy viscosity, p/p 

x, y, z physical coordinates 

X,Y,Z xtL,ytL,ztL 

X eddy viscosity parameter 

F normalized time scale 

|i viscosity 

5 boundary layer thickness 

p density 

(0 turbulence dissipation rate 

Subscripts 

/ freestream/inlet condition 

1 laminar flow 

o stagnation condition 

r reference condition 

t turbulent flow 

vv wall surface condition 

Computational Approach 

The numerical computational procedures used in this work 
were similar to the procedures in the NASA Glenn-HT 
RANS Code. The Code solved the compressible 3-D 
Reynolds-Averaged Navier-Stokes equations using a four 
stage Runge-Kutta scheme in conjunction with accelerating 
techniques (ref. 1). The Code had the capabilities of multi- 
grid and multi-block to accelerate the steady flow 
computations. In addition, a very low level of artificial 
dissipation was guaranteed by eigenvalue scaling. The Code 
used the k and co equations (ref. 2) to model the eddy 
viscosity. Instead of the existing turbulence model in NASA- 
Glenn-HT code, a one-equation turbulence model was used 
here to calculate the eddy viscosities to model the turbulence 
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terms in the Navier-Stokes equations. The current numerical 
computational approach is described in the following 
sections. 

Navier-Stokes Equations 

The compressible Navier-Stokes Equations in Cartesian 
coordinates without body force or external heat addition were 
used to describe the flow fields. The equations were written 
as 


dU dE dF dG n 
dt dx dy dz 


where 


(1) 
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U = 


p u 
pv 


p w 
E, 


t XX ^ (M 7 M- / )| 
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The Reynolds analogy 

k + k t =C p 


f vl+v± a 

V Pr / p h J 


E = 


p u 

pu 2 +p~T xx 
PUV-T X y 

puw — X x2 

(E, + p)u -in xx - vx X y - wx xz ~(k + k t )^- 


was used to relate the thermal conductivity and viscosity. Pr/ 
and Pr, were the laminar and the turbulent Prandtl numbers. 

The Navier-Stokes equations and their constituent 
equations were rewritten in a curvilinear coordinate system 
(ref. 1). These equations were then used in a finite-volume 
formulation (ref. 1) for numerical computations. 

One-Equation Turbulence Model 


F = 


pv 

P UV - X xy 

PV 2 + P — T yy 

PVW — X y- 

(Ef +p)v- ui xy - VT^ - wx yz ~(k + k t )^~ 


G= 


pw 

puw — X xz 

PVW — X y Z 

pw 2 + p- x zz 


{E, + p)w - ux xz - vx - wx zz - (k + k, )^~ 


dz 


The components of the shear stress were given by 


Assuming the original version (ref. 3) of the eddy viscosity 
transport equation was applicable for compressible flow, the 
author wrote the transport equation, in terms of a working 
variable v, as 


^-(pv)+^-(pMv)- 
ot ox 

t- J-(p' ; v)+^-(pwv 
dy oz 

_ 8 

p(v + v) dv 

d 

H 

p(v + v) dv 

dx 

a dx _ 

dy 

a dy 

d 

+ — 

dz 

p(v + v) dv 
a dz _ 

+ pC M (l - fti)Sv 

~ p{c w ifw - C bl f t2 / K 2 \v/df + p/ (1 Am 


The above equation omitted the contribution (ref. 3) of 
(grad v) 2 to the turbulence modeling and substituted with v 
for the diffusion coefficient (ref. 3). The eddy viscosity, 
V/ = |i/ p, was related to the working variable, v , through the 
equation 


V/ — y/vt 


NASA/TM— 2005-2 13830 


3 



where / v i = % 3 /(y 3 + C v i 3 ), / = v/v and v is laminar 
kinematic viscosity. 

The production term, S , is given by 

S = 5 + (v /K 2 d 2 )f v2 

S is the magnitude of the vorticity, d is the distance to the 
wall, and 

fv 2 =l-X/(l + X/vl) 


The destmction function, f w , is given by 


c 

fw — S 

v 


1 + C w3 6 

g 6 + C w 3 
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(2b) 


('-/-di-x/ft+xc,,)] 

Rez d 2 k 2 


where g = r + -r) and r = v /(5 k 2 z/ 2 ). 

The transition functions are 

.fa = C n g, exp(- C t2 [(o t 2 / A u 2 \d 2 + g, 2 d 2 )) 
ft 2 = Q3exp(-C M x 2 ) 

where g, = min(0. 1, At/ / co^Ax) . 

In the transition functions, co, is the vorticity at the surface 
and Ax is the grid spacing along the x direction. The velocity 
difference between a field point and the surface was Am. 

Spalart and Allmaras (ref. 3) used the following set of 
empirical constants: 

0,1 =0.1355,0,2 =0.622, ct = 2/3,k = 0.41 
C wl = C M / k 2 + (l + C b2 )/ a, C w2 = 0.3, C w3 = 2, 

C v i =7.1, Oi =1,02 =2,03 =1-1,04 =2. 

The turbulence transport equation, equation (2a) has the 
form similar to the compressible flow Navier-Stokes 
equations. The present author implemented equation (2a) in 
the Glenn-HT code to model the turbulence terms in the 3-D 
compressible Navier-Stokes equations. Similar to the 
manipulations of the Navier-Stoke equations in the Glenn 
HT-Code, equation (2a) was first non-dimensionalized with 
reference properties at the stagnation conditions: 

T r = T 0 ,p r = Po Ur = A,- = V 

and a length scale L. Equation (2a) becomes 


where p, v, d, and v are dimensionless quantities (with 
respect to the reference conditions and L ) and 


/ - 

- A 

2 

/ ^ - 


2 

/ - 

^ - \ 

dW 

dV 
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dU 

dW 

1 

dV 

dU 


~8Z) 

I 

[dZ 

~8X) 

I 

[dX 

~8Y) 


Equation (2b) showed that Re/ influenced the turbulence 
modeling. This was consistent with the dependence of C bl , 
C/,2, and C,i on the Reynolds number (ref. 3). The present 
author further assumed that the empirical constants, 
C/2, C/3, C/4, C w 2 , C v i,k, and a in equation (2b) had the same 

values as they were in reference 3. The values for the 
empirical constants, Cm, C/, 2 , and C,i were calibrated (ref. 4) 
for turbulent flat plate boundary layer flow. Computations, 
with C/,i = 0.1 or 0.16, C/, 2 = 0 or 0.622 and C,i = 0.1 did not 
predict a significant change in the U values. However, the 
computations, with C,]=I0 5 , predicted a reduction of the 
mean velocity at the near surface locations and predicted 
better the skin friction factors (C/= 3x10 ’’) at downstream 
stations. Although C,i had a small value, Re/xC,i could have 
a large value and it was consistent with the assembling of the 
turbulence transport equation (ref. 3). 

The eddy-viscosity equation was also used in the 
computations of the flows developed along concave wall 
surface without and with downstream flat plate wall surface. 
However, the following set of empirical constants, 

C M = 0.160, C/,2 = 0.622, o = 2/3, k = 0.41, 

C wi = Cm/k + (1 + C/, 2 )/o, C w2 = 0.3, C w 3 = 2. , 

C vl = 7.1, C n = 0.00001, C, 2 = 2„ C a = 1.1, C, 4 = 2, 

was assumed in the computations. For the computations of 
the boundary layer flows along the concave wall surfaces, 
C/3 = 0 was assumed for turbulence modeling (ref. 3). 
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The eddy viscosity modeling equation (2b) was rewritten 
in the curvilinear coordinate system as it was applied to the 
Navier-Stokes equations. The Navier-Stokes solver of the 
NASA Glenn-HT Code was updated to solve simultaneously 
the Navier-Stokes Equations and the eddy viscosity transport 
equation. 

Numerical Schemes 

Accurate and efficient numerical computational schemes 
to solve the Navier-Stokes equations were developed (ref. 1) 
for the NASA Glenn-HT code. The present author only used 
some of the numerical schemes in this work. The schemes 
used here are briefly described in below. 

A finite volume approach was used to discretize the 
governing equations. Fourth order cell center spatial 
discretization scheme was used in the present computations. 
The domain was divided into hexahedrons (cells). Flow 
variables were stored at the cell centers. On each cell face the 
convective and diffusive fluxes were calculated after 
computing the necessary quantities at the face center. Those 
quantities were obtained by averaging of adjacent cell-center 
values of the dependent variables. The system of discretized 
equations was advanced in time using an explicit four stage 
Runge-Kutta scheme (ref. 1) to reach steady state solutions 
of the flow variables. The local time step limit. At, was 
computed with both the convective, A t c and diffusive, A t d 
contributions as At = C 0 (A t c AtJ(At c + A t d )). The 
mathematical expressions for the convective and diffusive 
contributions could be found in reference 1. C 0 was a 
constant and was taken to be Courant-Friedrich-Levy (CFL) 
number. CFL = 1 was used in all the computations here. The 
maximum of the local time steps within the domain was used 
in the time-marching scheme to facilitate the steady state 
calculations. Artificial dissipation and residual smoothing 
(ref. 1) were also employed in the computations. The present 
author used one block and single grid in the computations. 
The multi-grid and multi-block options were not considered 
in this work. The present author used 2 different 
computational mesh configurations respectively in each flow 
computation and investigated the effect of the computational 
mesh configuration on the computations. 

Boundary Conditions 

Although three different flow fields were considered in 
this work, all the computational domains consisted of inlet 
and exit planes, wall surfaces and symmetric planes. Since a 
fourth order cell-center scheme was used for spatial 
discretization in the numerical schemes. The numerical 
computations required flow properties at two phantom cell 
centers immediately outside of the domain surfaces. These 
properties were specified with the following approaches. 

It was assumed that the cross flow velocity components 
were zero at the inlet phantom cell centers. Their static 
pressures and mass flow rates were set to the first interior 


cell centers’ static pressures and mass flow rates computed 
from previous time iteration. The density and total energy 
were then computed from the pressure and velocity with 
isentropic flow relationships. A uniform turbulent eddy 
viscosity parameter of the order of 0.1 was imposed at the 
inlet phantom cell enters. 

For the computations of the flows over the concave wall 
surfaces, flow symmetry was assumed at the far field domain 
surface. The flow properties at the phantom cell center next 
to the plane of symmetry were simply set to their properties 
at the interior cell centers adjacent to the boundary surface. 

A zero normal pressure gradient was assumed at the wall 
surface and the linearization approach (ref. 11) was used to 
set the velocity and the density at the phantom cell centers. 
The constant wall temperature condition was used in 
computations and the wall surface temperature was 4 percent 
larger than the room temperature. 

Flow symmetry was assumed at the domain surface along 
the Z direction and at the far field domain surface along the 
N direction. The flow properties at the Phantom cell centers 
next to these domain surfaces were set to their values at the 
interior cell centers adjacent to the boundary surfaces. 

A uniform static pressure was imposed at the exit plane. 
The other mean flow properties at the down stream phantom 
cell centers were extrapolated (ref. 11) from interior values 
and the exit static pressure. The exit pressure level for each 
flow field considered here is given in the following Results 
and Discussion section. It was further assumed that the 
turbulent eddy viscosity values at the downstream phantom 
cell centers were the same as the eddy viscosity value at the 
first interior cell next to the exit plane. 

Results and Discussions 

The numerical computations provided mean velocity 
components, mean temperatures, and eddy viscosities within 
3 different flows. The velocity, temperature, and eddy 
viscosity were studied for the predictions of the existing 
measurements of near wall surface velocity, the Reynolds 
shear stress, friction factor, and wall surface heat transfer 
rate. The following sections describe the computational 
predictions and their comparisons with measurements. 

Flat Plate Turbulent Boundary Layer Flow 

The computational procedures were used previously 
(ref. 4) to investigate the occurrence of a fully developed 
turbulent boundary layer velocity profile within a flow over a 
flat plate surface. The computations predicted approximately 
the law of the wall or the l/7 ,h velocity profile at a distance 
of 21 L ( L = 0.833 ft) from the plate leading edge. The 
computed velocity and eddy viscosity profiles at this stream- 
wise station were used here to study the calculation of the 
Reynolds shear stress profile within a fully developed 
turbulent boundary layer. 
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Figure 1. — Schematic of domain, 0-1-2-3-4-5-6-7. 


Some of the important parameters of the previous 
computations were summarized here. The computations were 
performed with one physical domain, figure 1. The domain 
had a physical dimension of 42, 0.7, and 0.7 along the X, 7, 
and Z directions. 


TABLE I.— MESH CONFIGURATIONS FOR FLAT PLATE 
BOUNDARY LAYER FLOW COMPUTATIONS. 


Mesh 

Domain 

Node 

Node 

ID 

Size 

Size 

Packing at 

1 

42, 0.7, 0.7 

120, 100, 6 

7=0, 0.7 

2 

42, 0.7, 0.7 

120, 100,6 

0 < 7< 0.05 


Two different computational meshes, table I, were used 
respectively to perform the computations. The meshes had 
the same node dimensions (120 in X direction, 100 in Y 
direction, and 6 in Z direction). The nodes were equally 
distributed along the X and Z directions but the meshes had 
different node distributions along the 7-direction. In mesh 1, 
the nodes along the 7 direction were packed at the near wall 
(7= 0) and free-stream locations (7= 0.7). In Mesh 2, the 
nodes were packed only at near wall (0 < 7< 0.05) locations. 
The free stream Reynolds number was 3xl0 5 /ft to ensure the 
occurrence of a turbulent boundary layer flow. With either 
mesh configuration, the computations predicted a skin 
friction factor, Cf « 3x10 at X = 21. The free-stream 
Reynolds number at X = 21 was 5.2xl0 6 . This Cf value 
agreed with the friction factor of a fully developed turbulent 
boundary layer flow. The profiles of mean velocity and eddy 
viscosity variations along the 7 direction at X = 21 were 
shown in figures 2a and b. The velocity profile indicated a 
boundary layer thickness, 5 of 0.35 L at this stream-wise 
station. This boundary layer thickness was the same as it was 
evaluated with the correlation of the boundary layer 
thickness and Reynolds number (ref. 12). The computed 
velocity profiles were also compared with the Law of the 
Wall and the l/7 th Law in figure 2b. 

In the present work, the author used the mean velocity and 
eddy viscosity profiles in figures 2a and b to calculate the 

Reynolds shear stress component, -mV at locations across 
the boundary layer thickness. The Reynolds shear stress was 
computed from 



0.0 50.0 100.0 150.0 

eddy viscosity, v t / V t0 


Figure 2a. — Boundary layer eddy viscosity 
profiles at x/L = 21. 



0.000 0.030 0.060 


Streamwise velocity, u/U f 
Figure 2b. — Boundary layer velocity 
profiles at x/L = 21. 



0.0000 0.0010 0.0020 
Reynolds stress, — u'v'jlJ j 

Figure 2c. — Boundary layer Reynolds shear 
stress profiles at x/L = 21. 
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-u'v' = v t (du/dy + dv/dx) (3) 

Figure 2c shows the computed Reynolds shear stress 
profiles in terms of the dimensionless parameter, 

2 

-u'v'/Uf , y/8 and mesh configurations. Experimental 
results (ref. 13) of the Reynolds shear stress within a 
turbulent boundary layer flow, with Re x = 4.2xl0 6 , were also 
shown in figure 2c. 

The results in figures 2a, b, and c show that computations, 
with two different meshes, predicted the same u, v ( , and 

-u'v'/u f at near wall surface locations (y/S < 0.05). The 
computations, using mesh 1 predicted rapid growth in eddy 
viscosities (fig. 2a) at 0.05 < yl 5 < 0.5 locations. The large 
eddy viscosity was associated with low mean velocities 
(fig. 2b). The velocity profile was close to the l/7 th law of the 
wall profile. Large Reynolds shear stresses (fig. 2c) were 
calculated according to the velocity and eddy viscosity 
profiles. The computed Reynolds shear stresses agreed well 
with their measurements (ref. 13). The computations, using 
mesh 2, predicted a low eddy viscosity (fig. 2a) and a high 
velocity (fig. 2b) within the boundary layer. The velocity 
profile at yl 5 >0.1 agreed well with the law of the wall 
profile (fig. 2b). The computed Reynolds shear stresses 
(fig. 2c) at y/8 > 0.05 locations decreased rapidly along the 
y direction and the shear stresses were smaller than the 
measurements. The results in figures 2a, b, and c showed that 
mesh configuration could affect the eddy viscosity, velocity, 
and Reynolds shear stress computations. 

Flow Over Concave Wall 

The author previously (refs. 5 and 6) used the 
computational methods to simulate existing experiments of 
free-stream turbulence and concave curvature effects on 
heated and transitional boundary layers (refs. 7 and 8). A 
sketch of the experimental flow filed and computational 
domain was shown in figure 3. The height, h, was 0.375 ft. 
The previous computational work showed the effect of 
computational mesh on the computational approach to 
predict the mean flow properties. The present work used a 
fine computational mesh to further explore the accurate 
predictions of near wall surface mean velocity, near surface 
temperature, and the Reynolds stress. 

In the present work, the computations were performed 
with a fine mesh which has 205 nodes in S direction, 
101 nodes in N direction and 57 nodes in Z direction. Along 
N direction, 33 nodes were packed in 0 < N < 0.01, the first 
node from the wall surface was at N = 0.0001 location and 
the node spacing was 0.0175 at N = 0.5. The inlet turbulence 
parameter, and the downstream static pressure, p d were 
0.075 and 0.96 respectively and these values were the same 
as they were in the previous work (ref. 6). Some of the 
present computational results were than compared with 


similar results from the previous work (ref. 6) with a 
different computational mesh configuration (mesh 4). 
Table II compares some of the important parameters between 
the two computational mesh configurations. 


TABLE II.— MESH CONFIGURATIONS FOR 
FLOW OVER CONCAVE WALL SURFACE. 


Mesh 

ID 

Packing 

5=0 

Packing 

5=12.5 

Packing 
N= 0 

Packing 
N= 0.5 

3 

0.059 

0.089 

0.0001 

0.0175 

4 

0.059 

0.089 

0.0001 

0.0300 

3 

Nodes 

205 in 5 

101 in N 

57 inZ 

4 

Nodes 

205 in 5 

61 in TV 

57 inZ 


With mesh 3, the computed values of skin friction factors 
converged when the computations proceeded to about 20,000 
steps of time iteration. As an example, the Cj values at 
3 different iteration steps, n = 19,800, 20,000, and 20,200 
and at locations along the S direction at Z = 0.044 are shown 
in figure 4a. The C/ values at four S stations from existing 
experiments (ref. 7) are also plotted in figure 4a to verify the 
accuracy of the computed C/ values. The computed C/ values 
jumped to large values at S = 3 station. Experiments (ref. 7) 
showed that flow transition occurred at this stream-wise 
station. Thus, the present computations with mesh 3 
predicted the transition occurrence and reduced Q values to 
their experimental values at greater downstream stations 
(S>4). 





Figure 3. — Schematic of flow over concave wall 
and computational domain, 0-3-4-5-6-7-2-1. 
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To investigate the mesh effect on the computations, the 
computed Cf values from the previous computations (ref. 6) 
with mesh 4 are shown in figure 4b. Mesh 4 had 205 nodes 
in S direction and 57 nodes in Z direction. There were only 
61 nodes in N direction and 23 nodes were packed in 0 > N > 
0.01 region. The second node was 0.0001 from the wall 
surface. The node spacing was stretched to 0.03 at far field 
boundary (N = 0.5). Figure 4b showed that the computed 
Cf values at 2 < S < 10 stations oscillated at n = 13,700, 
14,200, and 14,650 steps of time iteration. Mesh 4 had less 
nodes along the N direction and the computations did not 
predict a jump in the Cf value at S = 3 station. Thus, the 
computational approach required a relatively large numbers 
of nodes along the N direction to predict accurately the flow 
transition location and stable friction factors at locations 
along the stream-wise direction. 

The computed Q values (with mesh 3) at locations along 
S direction at n = 20,000 and at 3 span-wise Z stations 
(Z = 0.06, 0.044, and 0.088) are shown in figure 5. The 
Cf values peaked at S = 3 and 4 stations. The Cf values varied 
periodically along the Z direction at most S stations. Since Cf 
was computed from V s , the results in figure 5 indicated that 
V s at the near wall surface locations could vary periodically 
along Z direction. 

Figures 6a and b show the V s variations at N = 0.0025 and 
0 < Z< 0.1 locations at S = 3.12 and S = 5.35. The computed 
V s values with mesh 3 and at n = 20,000, are presented in 
figure 6a and the computed V„ with mesh 4 and at 
n = 14,650, are presented in figure 6b. These results were 
also compared with the measurements of the velocity 
components (ref. 7). Figure 6a shows that, at S = 3.12, the 
computed V s agrees well with the V s measurements and at 
S = 5.35, the computed V s varied periodically at the span- 
wise locations. The measurements did not clearly shows the 
span-wise periodic V s variation at S = 5.35 station. Figure 6a 
also shows good agreement among the computed and 
experimental V s values at Z = 0.044. This was consistent with 
the results of friction factor presented in figure 4a. 
Comparisons between the results in figures 6a and b shows 
that computational mesh configuration affected the accuracy 
of the V s computations. Computations with 205, 101, and 57 
mesh dimensions needed 20,000 time iteration steps to reach 
a converged solution. This simulated more closely the near 
wall V s span-wise variations at S = 3.12 and 5.35 stations. 

The computations used a constant temperature, 
TJT r = 1.04, as the temperature boundary condition at the 
wall surface. This wall temperature boundary condition was 
the same as the wall temperature of the experiment (ref. 7). 
The experiments applied a constant heat flux, Q w = 150 
W/m 2 to maintain the wall surface temperature condition. 

Assuming the one-dimensional heat conduction relation 
was valid between the wall surface and the first 
computational cell center locations, the surface heat transfer 
rate is given with 

Qw = -K (dT / 8N ) jv = o (4) 



Streamwise distance, S 


Figure 4a. — Friction factors computed with mesh 3. 



Streamwise distance, S 


Figure 4b. — Friction factors computed with mesh 4. 



Figure 5. — Friction factors at stream-wise stations, mesh 3 
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Figure 6a. — Spanwise ^variations at #=0.0025, 

S = 3.12, 5.35, n = 20,000 and mesh 3. 
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Figure 6b. — Spanwise ^variations at #=0.0025, 

S = 3.12, 5.35, n = 14,650 and mesh 4. 


Assuming that the temperature varied linearly at locations 
between the wall surface and the first cell center locations, 
equation (4) and the experimental Q w value were used to 
calculate the first cell center (at # = 0.00005) temperature. 
The resulting cell center temperature was used to verify the 
accuracy of the first cell center temperature from the Navier- 
Stokes calculations with mesh 3. Since the Navier-Stokes 
computations predicted closely the measured C/ values at 
n = 20,000 and at S = 3.2, 5.5, 7.7, and 9.9 stations (fig. 4a), 
the computed static temperatures at the first cell centers at 
those 4 S stations and at Z = 0.044 were investigated 





Streamwise distance, S 

Figure 7. — Near wall temperature and flow rate at 
Z = 0.044, n = 20,000 with mesh 3. 


according to equation (4) and the Q w value. The static 
temperatures at those cell centers obtained from the Navier- 
Stokes computations at n = 20,000 and from equation (4) 
were compared in figure 7. The first cell center temperatures, 
from those two different approaches, were about the same at 
the four S stations. The temperature within the flow field also 
depended on the mass flow rate conservation at stream-wise 
stations. The mass flow rates, 

0.5 

m = Jp V s dN (5) 

o 

within the S and # plane at S = 3.2, 5.5, 7.7, and 9.9 stations 
and at Z = 0.044 were computed and their values were also 
plotted in figure 7. The mass flow rate was well preserved at 
these S stations. The computations used the Reynolds 
analogy to model the turbulence contribution to the energy 
equation. Experiments (ref. 7) showed that the turbulent 
Prandtl number varied along the # direction at .S' = 5.3 and 
9.9 stations. Based on the Prandtl number profiles, an 
averaged Prandtl number = 0.75 was used in the present 
Navier-stokes computations. The effect of different Prandtl 
number on the computations was not investigated here. 

Experiments (ref. 7) also measured the Reynolds shear 
stress component, at locations along # direction at S = 5.3 
and 9.9 stations. The present author used these measurements 
to verify the accuracy of the Reynolds stress simulations in 
the present Navier-Stokes computation. The Reynolds shear 
stress component at each cell center was computed using 
equation (3) and velocity components at each cell center 
location. The computed Reynolds shear stresses, -u'v' at the 
cell centers along direction at n = 20,000, Z = 0.044 and at 
the S = 5.3 and 9.9 were plotted in figure 8a. These results 
were also compared, in the figure, with the measured shear 
stresses (ref. 7). The computed stresses were very small at 
S = 5.3 station and did not show the rapid decrease in the 
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Figure 8a. — Reynolds shear stress profiles at 
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0.35h 



Figure 9. — Schematic of flow over concave wall with 
extension wall and domain 0-3-4-5-6-7-2-1. 

Flow Over Concave Wall With Extension Flat Wall 


Distance from wall surface, N/8.5 
Figure 8b. — V s and V„ profiles along N direction, 

S = 9.9, Z = 0.044 and mesh 3. 

shear stress within a small distance from the wall surface. 
The computed shear stresses at locations along the 
N direction at S = 9.9 agreed well with their measurements. 
The measurements showed the effect of the direction of the 
normal velocity component, V,„ on the shear stress values. 
The downwash (V„ pointed toward the wall surface) reduced 
the shear stress value. The computed tangential and normal 
velocity components, V s and V„ at the cell center locations 
along the N direction at S = 9.9 are plotted in figure 8b. The 
computations predicted downwash ( V n < 0) at N! 8.5 < 0.012 
and upwash ( V„ > 0) at N/ 8.5 > 0.012 locations. The 
computed Reynolds shear stresses (fig. 8a) at the near surface 
locations reduced to the measured values with the downwash 
condition and then the computed shear stress increased, at 
locations away from the surface, to the measured values with 
upwash conditions. 


A sketch of the flow field and corresponding physical 
domain for the computations was shown in figure 9. The 
present computations used a computational mesh 
configuration (mesh 5) which had node dimensions of 405 in 
X direction, 91 in A direction, and 37 in Z direction. The inlet 
turbulence parameter, was 0.125 and the downstream 
static pressure, p d was 0.93 which were also used in the 
previous computations (ref. 6) with a different computational 
mesh (mesh 6). Table III listed a summary of the important 
parameters of the two computational meshes. A similar node 
spacing was assigned at the domain boundaries in both mesh 
configurations but large number of nodes was used in mesh 5 
to reduce the node spacing along the N direction. 


TABLE III.— NODE CONFIGURATIONS 
OF MESHES 5 AND 6. 


Mesh 

ID 

Node 

Size 

AN 

at 

N= 0 

AN 

at 

N= 0.5 

AS 

at 

5=0 

AS 

at 

5=24.5 

5 

405,91,37 

0.0001 

0.030 

0.030 

0.090 

6 

405,61,37 

0.0001 

0.030 

0.030 

0.090 
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The present computations (with mesh 5) showed that 
friction factors converged at 13,000 steps of time iteration. 
Then, computations were investigated for the accuracy to 
predict the skin friction factor and mean flow velocity, and 
the Reynolds shear stress modeling within the flow over the 
upstream concave wall surface. 

The computed skin friction factors, Cj at all 5 stations at 
Z = 0.01 and n = 13,000 are plotted in figure 10. Cf values 
obtained from previous computations (ref. 6) with mesh 6 are 
also plotted in this figure. The computations with mesh 5 or 
mesh 6 calculated almost the same Cf. The Cf value oscillated 
at downstream stations and the mean of C/was about 0.002. 
This mean value was close to the skin friction factor of a 
fully developed flat plate turbulent boundary layer flow. 
Experimental results (ref. 8) of the Cf values at 6 5 stations 
are shown in figure 10 for comparisons with the computed 
Cf values at corresponding 5 stations. The computed 
Cf values at near 5=10 stations agreed well with C, from 
experiment. The mean of the computed Cf at downstream 
stations was lower than the experimental Cf values at the 
downstream stations. 

Since the computations predicted the experimental 
Cf values at stations near 5=10, the author investigated 
further the computational results (with mesh 5) of the near 
wall surface mean velocities around the 5 = 10 station. 
Computed V s values at n = 13,000, N = 0.003, 0 < Z < 0.35 
and 5 = 7.7 and 9.9 are plotted in figure 11 and they are 
compared with the velocity measurements from the 
experiment (ref. 8). The experimental data was slightly 
shifted along the Z direction for better comparisons. 
Comparisons showed that computations predicted well the 
V s variation along the Z direction. Similar computational 
results from the previous computations (ref. 6) (with mesh 6) 
are shown in figure 12. The computed V s oscillated along 
Z direction especially at the 5 = 9.9 station. This was 
different from the computational results in figure 11. This 
showed the mesh effect on the V s computations. The 
computations with mesh 5 predicted more closely the 
experimental V s at near surface locations. 

The computed V s (with mesh 5) at cell centers along the 
N direction at 5 = 5.5, 7.7, and 9.9 are plotted in figure 13. 
Existing experiments (ref. 9) measured V s at locations along 
the N directions at 5 = 7.7 and 9.9 stations. The 
measurements are plotted in figure 13 to compare with the 
computed V s at 5 = 7.7 and 9.9 stations. Both the 
computations and experiments showed small differences in 
the V s values at N/ 8.5 > 0.002 locations. However, the 
computed V s values were larger than the V s values from the 
experiments. The computations also predicted a step increase 
within the V s profiles at 5 = 7.7, and 5 = 9.9 stations. The 
V s measurements at 5 = 7.7 and 9.9 did not show a step 
increase in V s values at large N locations. Figure 14 
compared the V s values at N < 0.01 locations and 5 = 7.7 and 
9.9 stations from the computations (with mesh 5 and mesh 6) 
and existing experiments (ref. 9). The computations 



Streamwise distance, S 

Figure 10. — Friction factors along S direction at 
Z= 0.01 with mesh 5 and mesh 6. 
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predicted well the experimental V s values and the 
computations and mesh 5 improved the predictions of the V s 
at N < 0.01 locations. The computations also showed that V s 
was reduced at S = 9.9 station. 

Equation (3) was also used to calculate the dimensionless 

Reynolds shear stress, u'v'/(dV s ) 2 , at locations along the 
N direction at S = 5.3, 7.7, and 9.9 stations. The profiles of 
the shear stress variations along the N direction at those 
S stations are shown in figure 15. Very small shear stress was 
calculated at S = 5.3 and 7.7 stations. The shear stress 
appeared at the S = 9.9 station. The maximum stress occurred 
at a location where V s started to increase sharply along the 
N direction (fig. 13). These V s and u'v'/(dV s ) 2 at large 

N locations agreed with experimental mean velocity and 
shear stress properties within a shear layer (ref. 14). 

As flow develops along the concave wall surface, vortices 
might appear within the flow field. The vortices could induce 
upwash (V„ > 0) and downwash (V n < 0) at the near wall 
surface locations (refs. 7 and 8). The downwash and upwash 
changed the V s profiles along the N direction at near wall 
surface locations. The vortices could lift off from the surface 
and induces shear layers at large N locations of downstream 
S stations. The computational results in figure 13 were 
consistent with the above argument. To investigate further 
the occurrence of vortices, the computed velocity 
components V s and V„ at N = 0.00005 (the first cell center 
from the wall surface), and two different Z locations at 
S = 53 stations are plotted in terms of the time iteration step 
in figures 16a and b. Figure 16a was from the computation 
with mesh 5 and figure 16b was from the computation with 
mesh 6. These results showed that both V s and V„ oscillated 
periodically at large time iteration steps (n « 10,000). V„ also 
changed its direction periodically with the time step. 
Downwash (V„ < 0) increased the V s value and upwash 
(V„> 0) decreased the V s value. 

To show some effect of downstream flat plate on the flow 
development along a concave wall surface, the computed 
velocity component V s and the Reynolds shear stress 

-mV/C/ f 2 at S = 9.9 for the cases without and with the 

downstream flat plate are compared in figures 17a and b. 
Measurements (refs. 8 and 9) of V s and the Reynolds shear 
stress for the case without downstream flat plate are also 
plotted in these figures. Figure 17a showed significant 
difference in the values and profiles of V s at locations along 
N direction. Different downstream static pressure boundary 
conditions were used respectively to perform the numerical 
computations for cases without and with downstream flat 
plate surface. The downstream pressure condition could 
change the V s variation. The downwash effect was significant 
for the case without the downstream flat plate and increased 
the Vs value at near surface locations. The V s profile also 
showed possible existence of a shear layer at 0.01 < 7V/8.5 < 
0.02 locations and the computations calculated an increase in 
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Figure 13. — V s profiles along N direction at 
S = 5.5, 7.7, 9.9, n = 13,000, and mesh 5. 
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Numbers of iteration, n 

Figure 16a. — Relationship of near wall V s and V„ 
at Z= 0.01, 0.175, S = 5.3, mesh 5. 



Numbers of iteration, n 

Figure 16b. — Relationship of near wall F, and V„ 
at Z= 0.08, 0.175, S= 5.3, mesh 6. 


the Reynolds shear stress, figure 17b, at these N locations. 
The free shear layer occurred at 0.03 < A/78.5 < 0.04 locations 
for the case with downstream flat plate. The downstream flat 
plate changed the flow development at locations near the 
concave wall surface and lifted the free shear layer further 
away from the wall surface. Figure 17b shows significant 
difference in the Reynolds shear stress predictions for the 
cases without and with the downstream flat plat. In 
turbulence modeling, the present study used different values 
for the empirical constant, C (3 , in the computations. C a = 1.1 
was used in the case without flat plate while C, 3 = 0 was used 
in the case with the flat plate at downstream. C, 3 = 0 was 
suggested (ref. 3) to model boundary layer flow 
development. The difference in the C /3 value might 
contribute to the difference in the computed flow properties 
in figures 17a and b. Figure 17b also showed the computed 



Figure 17a. — Comparisons of V s profiles without 
and with extension wall surface. 



Figure 17b. — Comparisons of Reynolds shear stresses 
without and with extension wall surface. 

Reynolds shear stress was small at small N locations. These 
could be due to the turbulence modeling assumption, 
v< = vf,i with /;,= % 3 /( 7 3 + C„i 3 ) and c vl = 7.1 x and f vl 
were small at the near surface locations which gave very 
small eddy viscosities at the near surface locations. 

This author also made a preliminary investigation of the 
inlet turbulence scale effect on the computed eddy viscosities 
at the near wall surface locations. The author used a low 
Reynolds number k-a> turbulence model (ref. 10) to compute 
the turbulence eddy viscosities at the near wall surface 
locations within the flow over the concave wall. The 
computations assumed 2-D flow exited over the concave wall 
surface. Only Mesh 4 was used in the computations. The 
computations were performed respectively with two different 
inlet turbulence scales, T\ = 0.00001 and 0.001 
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(non-dimensionalized with /?) and the same inlet turbulence 
intensity (ref. 6), T u = 0.006. The computed eddy viscosities 
at the near wall surface locations were compared with the 
eddy viscosities from the computations with One-Equation 
turbulence modeling. Figure 18a compared the turbulent 
eddy viscosities at N = 0.00015 and Z = 0.044 locations. 
With these two different turbulence models, the eddy 
viscosities had a level of 1.0 < v/v,/< 2.0 at all S locations. 
The eddy viscosities were close to each other at 3.5 < S < 8.5 
stations. The k-co modeling did not show significant inlet 
turbulence scale effect on the computed eddy viscosity. For 
clarity, only the results from the k-co with 7i = 0.001 are 
shown in figure 18a. Figure 18b shows the computed v„ k, 
and co variations, due to two different inlet turbulence length 
scales, along S direction at N = 0.00015 and Z = 0.044 
location. These results did not show a strong 7j effect on v„ 
k, and co values. These values also oscillated at 5.5 < S < 8.0 
stations. The oscillations might indicate the effect of flow 
transition. Flowever, the S locations were different from 
those from the experiment (ref. 6). This author only 
performed a preliminary study of k - co turbulence modeling 
for the eddy viscosity prediction. The results indicated that 
the inlet turbulence scale did not have a significant effect of 
the eddy viscosity at the near surface locations in the flow 
over the concave wall surface. 

Existing experiments (refs. 7 and 8) adjusted the lower 
wall surface geometry of curved duct flow to establish a 
boundary layer flow field along the concave wall surface. 
Measurements of the flow properties were taken within the 
boundary lay flow along the concave wall surface. The 
present study was directed toward the modeling of the flow 
properties at the locations near the concave wall surface, the 
effect of experimental lower convex wall contours on the 
flow development was not considered in the computations. 
From the experience of the present computations, the 
computational approach might require a huge mesh node size 
to model fully the experimental flow fields (refs. 7 and 8). 

Conclusions 

This paper presented the usage of a time marching 3-D 
compressible Navier-Stokes equation numerical solver and a 
one-equation eddy viscosity model to simulate the flow 
fields developed along concave wall surfaces with and 
without a downstream recovery flat wall. The Navier-Stokes 
numerical solver came from the NASA Glenn-FIT RANS 
code. The one-equation eddy viscosity model was derived 
from the Spalart and Allmaras model. Although 
computational approach, with the Spalart and Allmaras 
model and RANS solver, has been used to numerically 
simulate complex internal and external flows. The present 
author believes that the approach should be studied for 
simple flow field simulation before being widely used for 



Streamwise distance, S 


Figure 18a. — Comparisons of eddy viscosities with 
one-equation model and k - cb model. 



Streamwise distance, S 

Figure 18b. — Effect of turbulence length scales on k -(0 
model, Z = 0.044, N= 0.00015 and T u =0.006. 

complex flow analysis. Thus, this author performed the 
computational simulations of boundary layer flow 
development along concave wall surfaces. 

The author first described the governing equations and 
numerical schemes of the computations. The numerical 
computational schemes were validated for the predictions of 
Reynolds shear stress profile and velocity profile of a 
turbulent flat plate boundary layer flow. Then, the 
computational approach was studied for the prediction of 
surface friction factors, near wall surface velocity 
components, near wall surface temperatures and turbulent 
shear stresses within the boundary layer flows developed 
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along concave wall surfaces without and with downstream 
extension flat wall surface. The computational results of skin 
friction factor, near wall surface velocity component and 
temperature, and turbulent shear stress component were 
compared with their existing measurements and were 
discussed in terms of turbulence modeling, computational 
mesh configuration, and time iteration steps. The author 
concluded that: 

1. With a sufficient number of node points in a structured 
computational mesh, the computational approach 
predicted the occurrence of law of the wall mean 
velocity profile within a flow developed along a flat wall 
surface. The one-equation turbulence model could 
provide eddy viscosities to aid the computational 
methods to predict accurately a Reynolds shear stress 
component of a fully developed flat plate turbulent 
boundary layer flow. 

2. Given a mesh node packing condition along the normal 
direction to the concave wall surface, the computed skin 
friction factors showed the occurrence of flow transition 
at the upstream station. The computed skin friction 
factors and the near wall surface velocity components 
showed span-wise periodic variations as they were 
observed from the existing experiment. The 
computations predicted near wall surface temperatures 
which agreed well with their values derived from 
experimental wall heat transfer conditions. Based on the 
computational results of the mean velocity components 
and eddy viscosities, turbulence modeling of a Reynolds 
shear stress component agreed very well with shear 
stress measurements at downstream station. The 
computed shear stress also showed the effects of upwash 
and downwash on the shear stress. 

3. As a preliminary step to simulate a turbine blade flow 
with the present computational approach, it was used to 
compute a developing boundary layer along a concave 
wall surface with a downstream extension flat wall 
surface. The computational approach, with a fine 
computational mesh, predicted the experimental 
measurements of the near wall velocity components 
better. The computed near wall surface velocity 
components at the upstream stations showed the vortices 
effects on the upwash and downwash relationship among 
the velocity components. The computed velocity 
components at locations far away from the wall surface 
and at downstream stations showed the occurrence of 
shear layer due to the vortices lift off from the wall 
surface at the upstream stations. Based on the computed 
eddy viscosity and the velocity profile, the computed 
Reynolds shear stress was within the range of Reynolds 
stress measurements from existing shear layer 
experiment. 

4. With the aid of a low Reynolds number k- 0 ) turbulence 
modeling, the inlet turbulence length scale was found to 
have little effect on the computations of the eddy 


viscosities at locations near the concave wall surface. 
The eddy viscosities from the present One-Equation 
turbulence modeling were comparable with the eddy 
viscosities from the k-co modeling. Thus, the present 
One-equation turbulence model could be an alternate 
approach for turbulence modeling within flow near a 
solid wall surface. 

The proposed computational approach could include the 
effect of inlet turbulence level and concave wall curvature on 
the simulations of developing boundary layer flows. The 
inlet turbulence level and wall curvature play important roles 
in turbine blade flow development. Thus, the present 
approach may have the potential to improve the turbine blade 
flow analysis. Since single differential equation was used to 
model the turbulence in the Navier-Stokes computations, the 
present computational approach may reduce the 
computational resource requirement of the Navier-Stokes 
simulations of complex turbine blade flows. The turbulence 
model approach in this work also indicated the model 
dependence on the Reynolds number. The present author 
only performed preliminary studies of the Reynolds number 
dependency in the work. This turbulence modeling should be 
further verified with existing advanced turbulence modeling 
Finally, the present computational results showed the 
present computational approach could match the vortices 
effect on the near surface velocity in flow over concave wall 
surface. Vortices play important role in unsteady flow 
development. Thus, the present computational approach may 
lead to a foundation to develop a Navier-Stokes 
computational approach for the simulation of unsteady 
turbine blade flow field. 
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